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Imagine a scenario in which the dark energy forms via the condensation of dark matter at some 
low redshift. The Compton wavelength therefore changes from small to very large at the transition, 
unlike quintessence or metamorphosis. We study CMB, large scale structure, supernova and radio 
galaxy constraints on condensation by performing a 4 parameter likelihood analysis over the Hubble 
constant and the three parameters associated with Q, the condensate field: fig, Wf and zt (energy 
density and equation of state today, and redshift of transition). Condensation roughly interpolates 
between ACDM (for large zt) and sCDM (low zt) and provides a slightly better fit to the data 
than ACDM. We confirm that there is no degeneracy in the CMB between H and Zt and discuss 
the implications of late-time transitions for the Lyman-a forest. Finally we discuss the nonlinear 
phase of both condensation and metamorphosis, which is much more interesting than in standard 
quintessence models. 
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I. INTRODUCTION 

The observational evidence that the universe is acceler¬ 
ating today [1,2] changes the study of inflation from cos¬ 
mic archaeology to real-time ghost hunting and implies 
that our universe is dominated by energy with negative 
pressure or that Einstein’s General Relativity is incorrect 
on large scales. 

In this paper we discard the second possibility, but con¬ 
sider the idea that the dark energy comes from the con¬ 
densation of dark matter (which we take to be CDM). We 
show that the current data does not particularly favour 
the simplest model of condensation over the standard 
dark energy candidate: quintessence [3]. Nevertheless it 
has a very interesting nonlinear picture of the universe, as 
we describe in section V, and may serve as an archetype 
for other scenarios for dark energy. 

On a more general note, our analysis is perhaps the 
first detailed comparison with current data of a dark en¬ 
ergy model with a Compton wavelength, Ac ^ (U")“^/^, 
and speed of sound, = 6p/Sp, which change radically 
with time. By contrast, quintessence and metamorphosis 
[4] both involve scalar fields that are light at all times, 
though k-essence [5] by way of contrast, has a speed of 
sound which is low at decoupling, which may be observ¬ 
able in the CMB [6]. General issues with regard to the 
measurement of the dark energy speed of sound were dis¬ 
cussed in [7]. 

The idea that acceleration might be induced by a con¬ 
densate has been discussed by several authors. It has 
been suggested that a scalar condensate might arise in 
an effective field theory description of gravity [8]. Con¬ 
densates have also been invoked to explain inflation in 


earlier work [9,10]. These condensates would naturally 
form at very high energies, but from condensed matter 
theory we are also familiar with condensates forming at 
very low temperatures. 

An example is a super-conductive material where, be¬ 
low a threshold temperature, a phase transition takes 
place and the electrons (which are fermions) form Cooper 
pairs. Those Cooper pairs are collectively described by 
a boson field. In our cosmological scenario, when the 
temperature of the Universe reaches a suitable value, the 
cold dark matter, which could be associated to a fermion 
(e.g. the neutralino), might undergo a similar phase tran¬ 
sition where a condensate of fermion-fermion pairs would 
emerge. Such pairs would then be described by an effec¬ 
tive scalar field theory. Initially, the cold dark matter 
(CDM) scales as dust {wq = 0). If after the transition 
the final equation of state of the condensate Wf < —1/3, 
then the scalar held dynamics will drive the universe into 
an accelerating phase. 

The condensation mechanism has also been advocated 
recently in order to explain the late acceleration of the 
universe. In particular, in [19] it has been observed that 
a cosmological constant with the suitable energy scale 
(A ^ (10“^eV)‘^) can be explained by neutrino condensa¬ 
tion. Furthermore, in [20] it has been found that in super- 
symmetric, non-Abelian gauge theories with chiral/anti- 
chiral matter helds the condensation of those chiral helds 
in the low-temperature/strong-coupling regime yields a 
scalar held which has an inverse power-law potential and 
therefore is a possible quintessence held candidate. Sim¬ 
ilar ideas are explored in [21]. 

In an earlier paper we studied vacuum metamorphosis 
which also involves a late-time transition. We found that 
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condensation and metamorphosis & quintessence. Top: the 
Compton wavelength Ac = Bottom: The equa¬ 

tion of state w{z). In quintessence and metamorphosis w{z) 
has some slowly varying value during radiation and matter 
domination. In condensation, w{z) only vanishes on averag¬ 
ing over many oscillations, as for standard dark matter. After 
the transition at 2 = 2 t condensation becomes very similar to 
both metamorphosis and quintessence. For this example we 
chose Wf = —1. 


CMB, LSS and SN-Ia data prefer a transition around 
z ~ 1.5 — 2 [11]. One of the aims of this paper is to check 
whether the data simply prefer a late-time transition in 
the background cosmology (which controls the onset of 
acceleration) or are also fundamentally sensitive to the 
nature of the perturbations. To this end we compare 
metamorphosis and condensation. 

Within our phenomenological descriptions both of 
these have identical background dynamics. However, the 
perturbation physics is fundamentally different. In meta¬ 
morphosis the scalar field is very light (toq ^ lQ~^^eV) 
at all times*, both before and after the transition at Zt 
where w{z) changes. In condensation by contrast, the 


*In this paper we will denote the dark energy field by Q, 
whether it be quintessence, metamorphosis or condensation. 
Since the context should always be clear this should not lead 
to confusion. 


mass of the field is large before the transition [z > Zt). 
We choose m > IGel^ for z > Zt as appropriate for 
CDM, with a compton wavelength corresponding to nu¬ 
clear scales, Ac^lfm. Other choices, such as warm dark 
matter (m ~ IkeV) are certainly possible. After the 
transition (z < Zt) the scalar field becomes very light, 
with a Compton wavelength larger than lOOOMpc. We 
will show that the data are rather sensitive to these re¬ 
sulting differences in perturbation dynamics. 

There are at least four reasons why it is worth 
analysing a condensation mechanism for the origin of 
quintessence: Firstly, it can provide a link between dark 
energy and dark matter through a known physical phe¬ 
nomenon. This minimalist philosophy has been recently 
exploited in trying to obtain acceleration and clustering 
from the same source field [12] with a scale-dependent 
equation of state. It is also the basic driving idea behind 
the use of the Chaplygin gas which has equation of state 
p oc —p““, a > 0 and which interpolates between dark 
matter and dark energy [13-15], quite similar to the class 
of models which we study here. But it appears to have 
problems with age estimates of high-z objects [16]. 

Even more severe is that a recent analysis of fluctua¬ 
tions in Chaplygin-gas models finds that all models which 
differ significantly from effective ACDM are ruled out by 
the matter power spectrum [17] or the CMB [18]. Sand- 
vik et al find that if the sound velocity of the dark matter 
is not zero, its fluctuations are either unstable towards 
collapse (Cg < 0) or start to oscillate {cl > 0). Although 
the equation of state of the Chaplygin gas and our con¬ 
densation model are quite similar, it must be pointed out, 
that the condensation model is immune to this problem. 
The Compton wavelength of the dark energy becomes 
very large just at the transition (when the sound speed 
starts to differ from zero), and the small-wavelength fluc¬ 
tuations in this component are suppressed and cannot 
grow out of control. 

Secondly, condensation of CDM into a scalar field may 
give rise to scaling solutions characterised by a fixed frac¬ 
tion of the condensate relative to the CDM independent 
of redshift. Additionally the characteristic temperature 
of a condensation process is usually low, so such mecha¬ 
nisms naturally explains why a quintessence field would 
dominate only at low redshifts, although a condensation 
occurring in the future cannot be ruled out. In general a 
condensation mechanism for late-time acceleration would 
provide a partial solution to the coincidence problem. 

Thirdly, it is an example of a model where the effective 
mass of the dark energy field changes strongly. It allows 
us to study the influence of such a mechanism on CMB 
and LSS data, which probe the fluctuations as well as the 
background geometry. 

Finally, one interesting feature of condensation is its 
nonlinear origin. Hence, CDM condensation may give 
rise to significant modifications for the background cos¬ 
mology in regions of high density (such as clusters of 
galaxies) as compared to voids. Therefore, there exists 
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the possibility that condensation may be relevant to the 
core-cusp problem of CDM [23], an idea recently explored 
in [24], 

The aim of this paper is to analyse in detail how the 
condensation mechanism can be tested/falsified through 
its effects on CMB anisotropy, the matter power spec¬ 
trum and the supernova luminosity distance. We also 
compare this class of models with standard and non¬ 
standard quintessence models such as vacuum metamor¬ 
phosis [4]. The outline of the paper is the following: in 
the next section we describe the class of phenomeno¬ 
logical models which we have adopted for the conden¬ 
sation and how condensation of cold dark matter into 
a scalar field affects CMB anisotropies and the matter 
power spectrum. In section III we test the model pre¬ 
dictions with CMB, large scale structure and supernova 
data. In section IV we discuss how condensation differs 
from conventional quintessence and metamorphosis. We 
go on to examine how further constraints could be made 
with large scale structure analysis, angular-diameter dis¬ 
tance to radio sources measurements and detection of 
fluctuations in the dark energy. In section V we dis¬ 
cuss the non-linear evolutionary phase of the cosmologies. 
Section VI contains our conclusions and discussion. 


II. THE PHENOMENOLOGICAL MODELS 

In this paper we will compare and contrast three phe¬ 
nomenological models, particularly focusing on the first 
two: 

• Metamorphosis [4] - a light scalar field undergoes 
a transition in its equation of state from exactly 
pressure free dust to Wf < —0.3 at some redshift 
Zt- This is appropriate for approximately describ¬ 
ing many quintessence models [25]. 

• Condensation - A heavy scalar field with quadratic 
potential (e.g. CDM, dust on average) undergoes 
a transition in w{z) to Wf < —0.3 and a very light 
mass ~ I0“^^eV at Zt- Before the transition there 
are no light scalar fields or cosmological constant. 

• Evaporation - The inverse of condensation. At high 
redshifts z > Zt, there is no CDM. There is only a 
light scalar field and baryons. At Zt the scalar field 
condensate evaporates to yield the CDM, while 
w{z) becomes negative to drive acceleration at late 
times. 

We will now discuss in more detail these three cosmic 
scenarios. 


of the scalar field. Metamorphosis was studied in detail 
in [II] where it was found that the data prefers a late¬ 
time transition. The equation of state which we use to 
parametrise both metamorphosis and condensation is: 


w{z) = Wo + 


{wf - Wo) 

I -b exp((z - Zt)/A) 


( 2 . 1 ) 


where we assume that Wq = 0 and A = 30zt where A con¬ 
trols the width of the transition. This equation of state is 
then fed into a modified Boltzmann code which uses the 
reconstructed scalar field potential along the background 
trajectory of the field, Q(t) (see [11]). For an analysis of 
the effects of A on the CMB see [28]. For condensation 
this equation holds exactly for z < Zt but only on average 
for z > Zt- 

In terms of the scalar field dynamics, this class of mod¬ 
els is characterised by three free parameters: the redshift 
Zt at which the transition takes place, the fraction Qq of 
the total energy density today due to the scalar field and 
the final value w/ of the equation of state of the scalar 
field. 

We note that in fact this metamorphosis describes 
rather accurately a wide range of standard quintessence 
models [25], particularly the Albrecht-Skordis model [27]. 

t 

The main contrast between metamorphosis and con¬ 
densation is that in the metamorphosis case the Compton 
wavelength of the field Q is very large at all times, despite 
having w = 0 at early times. In condensation, as men¬ 
tioned earlier, before the transition (w) = 0 only, the 
Compton wavelength is very short until the transition, 
after which the field is very similar to the metamorpho¬ 
sis case; see fig. (1). 


B. Condensation 

The distinctive feature of condensate models is that the 
dark energy emerges from a condensation process which 
occurs when the universe reaches a certain critical tem¬ 
perature. This implies that above the redshift Zt at which 
the condensation occurs the dynamics of the background 
is determined by photons, baryons, CDM, and neutrinos, 
while below the redshift Zt cosmological dynamics is also 
determined by the dark energy field Q, its energy density 
PQ and its equation of state w{z). 

In general the details of the condensation will depend 
on the precise model used. We are not even guaranteed 
that such a condensate process is possible for realistic 


A. Metamorphosis 

The archetype of metamorphosis is the model proposed 
by Parker and Raval [4] in which non-perturbative quan¬ 
tum effects drive a transition in the equation of state 


^The parametrisation of w{z) in eq. (2.1) is interesting in 
terms of the statefinder pair {r, s} introduced in [29]. For a 
rapid transition, {r, s} = {1,0} for z > zt and ~ {1,0} for 
z < Zt it Wf ^ —1. The pair differ strongly from these values 
only near z = Zt where w{z) is large. 
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beyond-the-standard-model physics. Hence for simplic¬ 
ity and due to our ignorance of the detailed microphysics 
we shall assume that the condensation occurs instanta¬ 
neously on a constant energy hypersurface. Furthermore 
we impose that the total energy density is conserved 
through the transition. Since the transition is instan¬ 
taneous this implies that the total CDM-Q density per¬ 
turbation is conserved. 


PCDM ~ PCDM + Pq (2-2) 

^PCDM — ^Pq + ^PCDM (2-3) 

The fraction F of CDM energy density which is trans¬ 
ferred to the scalar field Q (pg = can be easily 

inferred from Eq. (2.2) 


F = 


fi 


Q 


^CDM 




1-f 




(2.4) 


where fig, flc dm are respectively the energy densities of 
the scalar field and CDM in units of the critical energy 
density today. The function is given by 


'i'wizt) 


3 r^* 

= fi Jo 


{! + *) 


™q(z) 


(2.5) 


In the following we shall consider a particular profile for 
the equation of state of the scalar field which coincides 
with that used for metamorphosis, i.e. equation (2.1). 

Using Eq. (2.4,2.5) it is straightforward to see how dark 
energy domination is addressed within condensation. To 
achieve ^q/^cdm = 0(1) implies that 

^ ( 2 . 6 ) 

If we assume that w{z) < 0 for z < Zj, condensation at 
high redshift has to produce a small fraction of scalar 
field energy density. Similarly, the more negative the 
equation of state of the scalar field is, the less efficient 
the condensation process has to be. 

In order to analyse the imprint of the condensation pro¬ 
cess on the CMB anisotropies and on the matter power 
spectra we need to solve the evolution equation for the 
scalar field fluctuations 5Q. The simplest initial con¬ 
ditions for the perturbations of the scalar field Q are 
5Q = 6Q = 0. However, since we are considering an 
instantaneous transition, this choice of initial conditions 
may not preserve adiabaticity In particular a discon¬ 
tinuity in the time-derivatives of the scalar field Q could 
generate a non-zero value of the intrinsic entropy pertur¬ 
bation of the scalar field, SqqI 


_ SQ SQ 


(2.7) 


This always vanishes on large scales [44] but could alter 
sub-Hubble dynamics if not treated carefully. 

In order to fix the initial conditions for the scalar field 
perturbations, we follow the method developed in [45] 
and we match the fluctuations of the scalar field Q along 
a constant energy density hypersurface, which appears 
the most appropriate for a condensation process. In par¬ 
ticular, the matching conditions for the perturbations are 
obtained by requiring the continuity of the induced met¬ 
ric and the extrinsic curvature. In the synchronous gauge 
and in a frame comoving with the CDM fluid, this im¬ 
plies that, for each Eourier mode dQk the quantity 0q 
defined as 

(pq(1 -I- wq))9q = —k^QSQk (2.8) 

must be continuous at the transition. Using the condi¬ 
tions (2.2,2.3) we have 0q = SQ = 0 across the transi¬ 
tion and SQ can be determined in terms of Q and 
by requiring Sqq = 0. From then on, the evolution 
of the perturbations is calculated using the same modi¬ 
fied Boltzmann code approach as for the metamorphosis 
model [11]. 


C. Evaporation 

Since we have considered the possibility that dark mat¬ 
ter condenses at some low temperature into dark energy, 
it seems natural to consider the inverse process. This 
would be the evaporation of a condensate to yield all or 
part of the dark matter of the universe. Less severe forms 
of evaporation have been studied within the context of 
Affleck-Dine baryogenesis [22] where a scalar condensate 
decays to give the required baryon asymmetry. 

In its extreme form (no dark matter before evapora¬ 
tion) and if the evaporation occurs suddenly and com¬ 
pletely at Zt then a clear disadvantage of evaporation is 
that the universe does not accelerate today since Dg = 0. 
It therefore makes sense to study evaporation models in 
which the transition has not completed by today and has 
a large width, A or not all the scalar field evaporates 
into dark matter. Can we put limits on such models? 
Let us consider the simplest model in which the there 
is no dark matter before the transition redshift Zt- For 
Wf < —0.5 and for z ^ Zt the universe is essentially a 
baryon-dominated flat universe since the condensate en¬ 
ergy density is negligible at high redshifts 

This places large constraints on the model since it is 
extremely difficult to reconcile CMB data, big bang nu¬ 
cleosynthesis (BBN), a flat universe and SN-Ia data (the 
incompatibility is 3 ct or more [58] in a A baryon domi¬ 
nated model). In addition, the large baryon content of 


* We assume the primordial spectrum of fluctuations is purely ^Though note that if w = 0 and Uq is not negligible at 
adiabatic. decoupling this will not be true. 
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the universe tends to cause very prominent oscillations 
in the matter power spectrum that are not strongly ob¬ 
served in current data. This latter test is less straight¬ 
forward to apply to evaporation since the nascent dark 
matter might tend to smooth out the oscillations if it is 
warm enough to have a sufficiently large free-streaming 
length. 

A very interesting point is that evaporation would lead 
to a strongly time-dependent bias around the transi¬ 
tion redshift, an effect which naively links to the work 
of Magliochetti et al [50] who found a strong redshift- 
dependence of the bias. 

From the results of Griffiths et al [58] we can conser¬ 
vatively infer that, for Wf ^ —1, a low zt is unacceptable 
for evaporation unless wq is close to 0. In fact, since the 
secondary acoustic peaks are sensitive to the baryon and 
dark matter content, it is likely that best fit to the CMB 
would occur for zt > 1000. 

Such values would be disfavoured from SN-Ia observa¬ 
tions unless the transition were extremely wide so that 
there still remained a significant component with neg¬ 
ative equation of state at low redshifts. As such, only 
’’glacier evaporation” models would seem to be observa- 
tionally viable. However, the evaporation process of the 
condensate in such a slow manner is physically realistic 
and evaporation has some interesting nonlinear implica¬ 
tions which we discuss in section V. We will not, however, 
explicitly compare evaporation models with current data, 
leaving that to future work. 

III. COMPARISON WITH OBSERVATIONS 
A. Data and analysis method 

In order to constrain/falsify our phenomenological 
models of condensation, we compare its predictions with 
a number of observations. In particular we consider three 
different kinds of datasets: CMB anisotropy, large scale 
structure (LSS) and type-la supernova (SN-Ia) data. For 
the CMB anisotropies, we use the decorrelated COBE- 
DMR data [33] and the recent data from the DASI [34], 
MAXIMA [35], BOOMERanC [36] and CBI [37] exper¬ 
iments. The total number of data points is ncMB = 58, 
with I ranging from 2 to 1235 for the first four sets, and 
the mosaic data from CBI extends this to 1900. In the 
analysis, we take into account calibration uncertainties 
for DASI, MAXIMA, BOOMERanC and CBI. 

As far as LSS data are concerned, we use the matter 
power spectrum estimated from the following surveys: 
2dF 100k redshift survey [38], IRAS PSCz 0.6 Jy [39], 
and Abell/ACO cluster survey [40]. In order to avoid 
possible effects due to non-linear contaminations, we fix 
a cut-off fcmax = 0.2/i/Mpc. The total number of data 
points is riLSS = 48. Einally, we use the redshift-binned 
SN-Ia data [42] from HZT and SCP (nsN = 7). Ad¬ 
ditionally we also studied the angular-diameter distance 




FIG. 2. A compilation of our data for the CMB (top) and 
LSS (bottom) used iu our analysis. The PSCz and 2df are 
galaxy surveys while the Abell/ACO points are for cluster 
data. The supernovae data are as described in [11]. 


to radio sources [41]. This method is potentially able 
to probe a large range of redshifts, current data already 
covers 0.1 < z < 4. All of our models as well as ACDM 
are perfectly consistent with current data, which pro¬ 
vides a valuable test. Unfortunately we will have to wait 
for larger samples before we can constrain our additional 
parameters significantly. 

We performed a likelihood analysis by varying 
the three parameters directly related to the scalar 
field {VlQ,Wf,zt) and the Hubble parameter, = 
lOOh-okm,s“^,Mpc“^. For the other parameters we 
chose the following values: = 0.02, Ug = 1, 

n* = 0, T = 0, Htot = 1) where r is the reionisation 
optical depth and are the spectral indices for the 
scalar and tensor perturbations respectively. We also in¬ 
clude the “standard” effective number number of neu¬ 
trinos (3.04) and we fix the cold dark matter density as 
f^CDM = 1 — <^b/(^o) ~ ^Q- la our analysis we set the 
amplitude of the tensor perturbations to zero. 

In order to compute the likelihood functions for 
the various parameters we follow the usual proce¬ 
dure. We compute the power spectra Cg of the CMB 
anisotropies and the CDM matter power spectra P{k) 
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for each set of parameters over the four-dimensional space 
{nQ,Zt,Wf,Ho). We then evaluate the ^it each grid 
point. The CMB anisotropy spectra and the matter 
power spectra are related through their respective ampli¬ 
tude normalisations. We found that the CMB data fixes 
quite well the overall amplitude for each model. As far as 
the LSS data are concerned we allowed a bias 1/5 < 5 < 5 
for 2dF and PSCz and 1/9 < 5 < 9 for Abell/ACO since 
clusters are expected to be more biased than galaxies. 
We choose these large bias values since the biasing mech¬ 
anism is still not very well understood especially in the 
context of non-standard dark matter models, and since 
(as argued in section V) the possibility of strongly non¬ 
linear physics enters on cluster scales. 

We then finally determine the one-dimensional and 
two-dimensional likelihoods by integrating over the other 
parameters. 


B. Results 

Our main results are shown in figures (3) and (4). Both 
models prefer a energy density Oq close to 0.7 (0.7 for 
metamorphosis, 0.75 for condensation). Both models also 
prefer a final value of the equation of state Wf < —0.8. 
A big difference is seen in the likelihood curve for zt- 
Whereas metamorphosis favours Zj ^ 1.5—2.5, in conden¬ 
sation a late transition {zt < 1.5) is strongly disfavoured. 

Overall the condensate model and the metamorphosis 
models both have slightly better best-fit than the best 
ACDM model, although this is to be expected since we 
have extra parameters and is not highly significant. The 
condensation best fit = 79.3, with parameters Zt = 4, 
Wf = —0.95 and Qq = 0.75, while the ACDM best fit 
X^ = 84.9 (for Da = 0.73). The best fit for metamorpho¬ 
sis has = 78.8, with parameters Zt = 1.5, Wf = —1 
and Dq = 0.73. 

The following sections discuss in detail the CMB and 
LSS results. It should be noted that the supernova data 
only probes the background geometry, which is the same 
for both metamorphosis and condensation. Since it is 
also independent of the Hubble constant, the results of 
[11] are unchanged, except that we have updated the 
brightness of the supernovae at z ^ 1.7 to account for 
lensing [43]. Namely, there is a slight preference for a 
low-z transition, but the relatively large error bars and 
lack of points at z > 1 do not allow strong conclusions. 
Section IV C discusses more general aspects of using dis¬ 
tance measurements to probe the background geometry 
of the universe. 


CONDENSATION 



FIG. 3. The Id-likelihood curves ( Top: Condensation, 
bottom: Metamorphosis) for the scalar field parameters 
(fig, zt,Wf) obtained combining all the CMB, LSS and SN-Ia 
data. The dashed curves are the likelihoods obtained by in¬ 
troducing a Gaussian prior on Hq {ho = 0.72 ± 0.08). 


IV. COMPARISON OF CONDENSATION AND 
METAMORPHOSIS 

A. CMB 

The likelihood plots for condensation and metamor¬ 
phosis are very similar when considering Dg and w /. The 
major difference comes with the redshift of the transition 
Zt- Although the background dynamics are identical, 
Metamorphosis favours a transition of z* ~ 2, whereas 
condensation strongly disfavours Zt < 1.5. 

To explain this we need to take a closer look at the 
respective CMB angular power spectra, see fig. 7. Meta¬ 
morphosis, like quintessence, has an almost trivial effect 
on the CMB relative to ACDM, simply through the ISW 
effect which changes the Sachs-Wolfe plateau and hence 
the relative height of the peaks through the CMB nor¬ 
malisation on large scales. The overall form of the CMB 
power spectrum of metamorphosis is not changed. 

The effect of Zt on the acoustic peaks of condensa¬ 
tion is more subtle. Although the total density of all 
species is unchanged, the phase-transition means that the 
change in the CDM density is not continuous. Changing 
Zt changes the CDM density at decoupling and the red- 
shift of matter-radiation equality, even though its value, 
ilcDM — 1 — Dg, is unchanged today. This alters the 
gravitational potential wells, although the effect is only 
significant for small Zj, where the amount of CDM that is 
condensing out is at a maximum. Decreasing the amount 
of CDM changes the time of matter-radiation equality, 
which boosts the amplitude of all oscillations through an 
ISW effect. The third peak is moved to slightly higher £ 
and thereby damped more, making it less high than the 
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CONDENSATION 



FIG. 4. The 2d-likelihood contours (Top: Condensation, 
bottom; Metamorphosis) for the scalar field parameters 
(yiQ,zt, Wf) obtained combining all the CMB, LSS and SN-Ia 
data. 


second peak [47-49]. Overall, the condensation CMB 
power spectrum effectively interpolates between the two 
extreme cases of sCDM (A = 0,11 = 1) and ACDM. 

For interest’s sake we also plot in hg (8) the predic¬ 
tions for the temperature-polarisation (TE) cross-power 
spectra for condensation and metamorphosis and com¬ 
pare them with the data from the first year of WMAP 
data [?]. Currently the data cannot distinguish between 
the models, especially as the curves are sensitive to reion¬ 
isation history which will in turn depend on zt etc... 

B. LSS 

The large scale structure data also strongly disfavours 
any value of the transition redshift zt < 2 for condensa¬ 
tion. In contrast, the case of metamorphosis, Zt seems to 
be almost independent of the LSS data, as shown in figs. 
(5,6). Why is this? 

Comparing the transfer functions, T{k), for the meta¬ 
morphosis and condensation models with ACDM and 
sCDM (where the universe is flat and CDM dominated 
without cosmological constant so Dcdm = 1 — Db) is in¬ 
formative. When normalised to give the correct COBE 
normalisation, sCDM gives a slightly smaller amount of 
power on large scales, but gives more power on small 
scales. 

In condensation the Compton wavelength is small be¬ 
fore the transition, and only becomes big after the tran¬ 
sition. Clustering therefore does take place on all scales 
for z > Zt- For z < zt, the growth of clustering on small 
scales ceases due to the large Ac but this also occurs on 
all scales since the universe starts to accelerate. The 
closer that transition is to today, the more clustering on 
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FIG. 5. Condensation: the Id-likelihood functions for 
the scalar field parameters (Dq, 2t, w/). Left-column: CMB. 
Middle column: LSS. Right column: SN-Ia. The thick solid 
orange curves are the likelihoods obtained by introducing a 
Gaussian prior on Hq {ho = 0.72 ± 0.08). 


CMB LSS SNla 



FIG. 6. Metamorphosis: the Id-likelihood functions for 
the scalar field parameters {i^Q, zt,Wf). Left-column: CMB. 
Middle column: LSS. Right column: SN-Ia. The thick solid 
orange curves are the likelihoods obtained by introducing a 
Gaussian prior on Hq {ho = 0.72 ± 0.08). 
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Condensation 


Condensation 




FIG. 7. The Ce spectra for condensation (top) and meta¬ 
morphosis (bottom) as a function of zt. For comparison in 
the condensation case sCDM {fl = 1,A = 0) is the lowest 
curve at £ = 200, ACDM the highest. The other parameters 
are fixed to be fig — Ha ~ 0.7, ho — 0.63 and Wf = —0.95. 


small scales that can occur, resulting in a larger mass 
dispersion, erg. Again, the power spectrum interpolates 
between sCDM {zt 0) and ACDM {zt ^ 1). Metamor¬ 
phosis, in contrast, is a light scalar held before and after 
the transition, with little change in the Compton wave¬ 
length and so the large scale structure is only weakly 
dependent on the transition redshift, as discussed in [11]. 

But why does the LSS data disfavour a low zt for con¬ 
densation? We have not computed as for all these models 
(we do so below for specihe parameter values), so the rea¬ 
son cannot be that. Similarly, our bias upper limits of 5 
and 9 for galaxy and clusters respectively are unlikely to 
have a real impact in excluding low-zj models. 

The key lies in the change to the redshift of the matter- 
radiation equality, Zgq, which changes signihcantly with 
Zt for Zt < 5 and provides the dominant effect on the 
shape of the power spectrum. Comparing the LSS data 
(fig. Ill A) and the condensation T{k) (fig. IV B) we see 
that the sCDM power spectrum turns over around k ~ 
0.1hMpc~^ while for zt = 0.5 the turnover starts around 
0.05hMpc~^. In comparison, the data do not show any 
definite sign of a turnover, even out to k = 0.01hMpc~^. 



Metamorphosis 



FIG. 8. The Cj^ temperature-polarisation spectra for con¬ 
densation (top) and metamorphosis (bottom) as a function of 
Zt vs the WMAP data. For comparison in the condensation 
case sGDM (12 = 1, A = 0) and AGDM curves are included. 


Hence both the CMB and the LSS disfavour condensation 
with Zt < 1.5. The same conclusion is likely if one does 
include as constraints. 

The overall normalisations of the power spectra are not 
trivial, since they are linked to the COBE normalisation 
of the CMB power spectrum. We can read it off the high- 
k region of fig. 9. While condensation lies between ACDM 
and sCDM as expected, we have to remember that the 
ISW effect of metamorphosis changes the normalisation 
strongly and leads to a low value of as [28]. 

The differences in perturbation evolution mentioned 
earlier can also be seen in fig. 10, which shows as (rel¬ 
ative to the one of a ACDM universe) as a function of 
redshift for zt = 1.5 and Wf = —0.95. Standard CDM 
starts to produce more structure, especially at low red- 
shifts where a cosmological constant would start domi¬ 
nating. Condensation behaves for z > zt like sCDM, but 
its behaviour changes at Zt roughly when the universe 
begins to accelerate, thereby forcing the linear perturba¬ 
tions to stop growing. 

The redshift dependence of as and the bias factor can 
therefore be used to distinguish these models. Maglio- 
chettiet al [50] found a strong signal of a redshift- 
dependent as for z > 2 and redshift dependent bias. The 
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CONDENSATION 



- CDM 

- ACDM 

. z=0.5 

- z = 1.5 


-z=0.5 

- z=1.0 

- z=1.5 

- z=3.0 


FIG. 9. The COBE normalised matter power spectra for 
condensation (top) and metamorphosis (bottom) as a fnnc- 
tion of zt, with sCDM and ACDM to compare with in 
the condensation case.The other parameters are fixed to be 
~ ~ 0.7, ho = 0.63 and Wf = —0.95. 



FIG. 10. The ratio of ag (mass dispersion smoothed over a 
radins of 8h~^ Mpc) as a function of redshift for the different 
models, as compared to ACDM. Da = 0.7 for the ACDM 
model, and 0 for sCDM. Dq = 0.7, Wf = —0.95 and zt = 1.5 
for the condensation and metamorphosis models. The effect 
of the growth of small-scale perturbations due to the small 
Compton wavelength in the condensation case, for z > zt is 
clear. 


redshift dependent bias b(z) is expected in our models 
since the baryons will take some time to respond to the 
sudden change in the dark matter characteristics. We 
leave this interesting issue to future work. 

As a final point we note the strong differences be¬ 
tween the CMB Id likelihood curves in figs. (5) and 
(6) with and without the HST prior on the Hubble con¬ 
stant. In particular, Wf changes drastically reflecting the 
well-known CMB degeneracy between H and Wf. 


C. Distance measurements 


One important tool to constrain the equation is state is 
the determination of distances as a function of redshift, 
e.g. with supernovae. The basic building block is the 
comoving distance (for a flat universe) 

dc{z) = c(77o - ri{z)) = — f 7 ^ (4.1) 

ao Jo H{z) 


where r]{z) is the conformal time and 


H{z) = Ho 


1/2 




(4.2) 


i indexes the constituents of the energy density in the 
universe, eg. radiation (w = 1/3), matter {w = 0), cur¬ 
vature {w = —1/3) and quintessence type contributions, 
where w can vary as a function of z like in our mod¬ 
els. The apparent distance to a bright object such as a 


supernova is then determined by the luminosity distance 
dp = {l + z)dc, while the distance to an object of a given 
physical length (like radio sources, see [41]) is described 
by the angular diameter distance dA = dc/(l + z) **. 
It is also possible to use the CMB in a similar way: the 
angle subtended by the sound horizon at recombination 
gives us a fixed scale, and we can use it to determine the 
geometry between and today. Since the horizon size 
is a comoving quantity and not a physical quantity, it 
is dc{zrec) which enters in this case. (See e.g. [59] and 
references therein for a detailed discussion.) 

While it seems at a first glance that the distance mea¬ 
sures depend on the value of the Hubble constant Hq, 
this is not so: all data available compares two lengths 
and thereby divides out any dependence on Hq. In the 
CMB case, the size of the sound horizon at last scatter¬ 
ing depends linearly on I/Hq just like dc; in measuring 
the luminosity distance to SN-Ia, the normalisation of 
the distance ladder is treated as a nuisance parameter 
and effectively determined in units of l//i, as is the case 
for the physical reference length when working with the 
angular size of radio sources. 

It is not possible to determine several quantities from 
only one data point (e.g. the location of the CMB 


**dL and (Ia are, of course, equivalent and related via the 
reciprocity relation. 
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FIG. 11. The left figure shows the degeneracies for mea¬ 
suring Q.m versus Q.a with supernovae at 2 = 1 (dashed lines) 
and the position of the CMB peaks (solid lines). The degen¬ 
eracies are clearly orthogonal and an accurate measurement 
of both variables is possible. The right hgure shows the same 
for w and Q.q and measurements at 2 = 1 (dotted), 3 (dashed) 
as well as the CMB peak position (solid). The likelihood is 
very flat in the w direction, making it difficult to pinpoint 
it precisely, especially for a very negative equation of state. 
The models shown here have a constant equation of state [zt 
is very large). 



Oa 


FIG. 12. The left hgure shows the degeneracies for a model 
with Wf = —1 and observations at 2 = 0.5 (dotted), 2 = 3 
(dashed) and the CMB (solid). The right hgure shows the 
same for fig =0.5. A priori it is necessary to bracket the 
phase transition. As the right plot shows, high -2 distance 
data can help hnding the transition in certain cases, but does 
usually not improve much on the CMB data (which is ehec- 
tively at 2 « 1100). Also, transitions at redshifts higher than 
5 will be very hard to detect unambiguously using distance 
measurements. 


peaks) Measuring distances over a range of 2 values as 
done for supernovae and radio galaxies, improves the sit¬ 
uation somewhat, as does combining different data sets. 
Indeed, a low zt can be strongly constrained by data on 
both sides of the transition (see figure 12). A high zt'>5 
may be impossible to detect given the current limits on 
w, since the dark energy becomes irrelevant very quickly. 

Figure 11 also reinforces what the full likelihood figures 
in section III have already shown: fig can be determined 
to a much better accuracy than w, since distance mea¬ 
surements depend only quite weakly on its value. 

It must be emphasised at this point that especially 
the CMB degeneracies depend strongly on the quantities 
which are being kept constant. If we fix instead of 

the Hubble constant (and therefore vary h as a function 
of flm), the curves look substantially different (see e.g. 
[57]). Nonetheless, the overall conclusions remain the 
same. 


D. Change in the comoving number density of the 
Lyman-o absorbers 

The number of Lyman a absorption lines per unit 
redshift, dn/dz, has been of interest to observers for a 
long time. It is sensitive both to the background evolu¬ 
tion and to details of the linear and nonlinear density 
perturbation evolution [51]. This makes it a possible 


^^Of course the full CMB power spectrum contains many 
more observables, thereby breaking some of the degeneracies. 


probe of quintessence type models. Unfortunately, the 
non-gravitational effects make this instrument very diffi¬ 
cult to analyse and interpret reliably and would require 
high-resolution hydrodynamic N-body simulations. Nev¬ 
ertheless we illustrate the underlying point with a simple 
model. 

Observers use a comoving volume element with a fixed 
physical radius Vp and comoving length dvc [52]. As Vc = 
cry, we find that (irc/(i 2 = 1 / 77 ( 2 ), and since Cc = (l-|- 2 )rp 
we find that 


dV = drl oc (1 -b zf/H{z)drldz. (4.3) 

If we naively assume that the comoving number density 
of Lyman a absorbers is constant, so that they directly 
probe the volume as a function of redshift, we set dn{z) = 
PodV ( 2 ) and obtain 


dn 

dz 


PoHo 


77(2) 


(4.4) 


Observers use the parametrisation dn/dz oc (1 -b z)z. A 
pure CDM universe would therefore have 7 = 1/2, while 
a pure A universe had a 7 = 2 . As observers find 7 > 2 
even at redshifts larger than 1 [53] we have to conclude 
that either the average equation of state is extremely 
negative at these redshifts, or that the non-gravitational 
evolution effects dominate. 

If it became possible to study the cosmological evolu¬ 
tion only, then the density of Ly-a absorbers would be an 
extremely useful tool since the change of physical density 
depends on H{z) and not its integral (as is the case for 
distance measurements and the position of the first peak 
in the CMB, see section IV C). Unfortunately the actual 
change of the number density is strongly influenced by 
evolutionary effects. At the moment one can speculate 
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FIG. 13. An illustration of the break in the exponent 7 de¬ 
rived from (4.4). While a ACDM universe shows a smooth 
transition from a steep slope at low redshift to a slope of 
7 = 1/2 when CDM starts to dominate (dashed line), a con¬ 
densation or metamorphosis model exhibits a break at the 
transition redshift (here = 1, solid line). Although Lyman 
a data is strongly influenced by non-gravitational evolution 
effects, we would expect to see some kind of jump in d'y/dz if 
a transition takes place at a low enough redshift. 


whether it is possible to draw conclusion from a differen¬ 
tial approach: figure 13 shows dn/dz for a ACDM model 
and for a universe with a transition at Zt = 1. The sec¬ 
ond case shows a break in the comoving number density, 
to a slope with a lower 7 . This change in the shape of 
dn/dz might be detectable even though the actual form 
of the curve is strongly modified by additional physics. 

As a final comment we note that there is also evidence 
for a change in the Lyman-a optical depth at a redshift 
z ~ 3.2 [54]. 

E. Detection of fluctuations in the dark energy 

Even though scalar field dynamics can lead to a back¬ 
ground evolution very close to the one of a ACDM cos¬ 
mology, there is a basic difference: A cosmological con¬ 
stant is a smooth quantity, while the scalar field can have 
fluctuations [11,55,28]. In the models considered here, 
these fluctuations can introduce modifications of the Cy 
spectrum of the order of a percent or more, even in the 
worst case of a large Compton wavelength at all times. 
Especially in the condensation scenario, it is impossible 
to avoid those fluctuations, since they are present in the 
dark matter. The matching conditions will transfer them 
into the dark energy component. An example is provided 
by the two curves for metamorphosis and condensation 
with zt = 0.5 in fig. 7. Both have the same effective 


equation of state, but the resulting CMB power spectra 
are very different. 

An experiment like Planck, which is basically cosmic 
variance limited at low to average £ {ae « l/v^2F+T), 
could detect (or rule out) these fluctuations for some of 
the more prominent cases (like a Zt <C Zlss or a small 
Ac). Even though this effect will be very difficult to dis¬ 
entangle from the influence of other cosmic parameters, 
it would provide an important insight into the nature of 
the dark energy. 

V. THE NONLINEAR EVOLUTIONARY PHASE 

Both the condensation and metamorphosis models 
have an extremely interesting nonlinear portrait. Since 
the exact time and location of the the two transitions is 
sensitive to local conditions of density and temperature, 
areas of nonlinearity can significantly delay the onset of 
the transition. 

Let us consider metamorphosis first. In metamorphosis 
the transition occurs when the Ricci scalar grows to be 
of order —m^/(^ — 1/6) [11]. Now R = —kT where T = 
p—3p ~ p for CDM. Hence, in high-density regions which 
will separate from the Hubble flow and recollapse to form 
bound structures like galaxy clusters, the Ricci scalar 
is more negative than in the low-density inter-galactic 
medium (IGM). 

This implies that the IGM will be the first place where 
the metamorphosis transition occurs. Inside rich clusters, 
it may even occur that the local nonlinear growth of den¬ 
sity perturbation overcomes the average FLRW growth 
of R towards zero, and causes it locally to begin decreas¬ 
ing again. In these regions the metamorphosis transition 
might never take place. 

The first implication of this is radical - the perturbed 
FLRW formalism for describing the universe rapidly be¬ 
comes invalid as the universe now consists of a fluid with 
an inhomogeneous background equation of state - nega¬ 
tive in voids and positive in rich clusters. 

This will lead rapidly to a strongly scale dependent bias 
since density perturbations in voids and in clusters will 
evolve very differently. On scales smaller than the inho¬ 
mogeneity in the equation of state, the perturbed FLRW 
approximation is good. Inside clusters perturbations will 
grow as they would in a flat, unaccelerating GDM model 
while in voids the perturbations will feel the expansion. 

This has important implications since the issue of clus¬ 
ter abundance versus redshift suddenly becomes an ex¬ 
tremely complex problem. It is no longer clear that clus¬ 
ters will be assembled at low-redshifts as happens in open 
or A-dominated universes due to the freeze-out in pertur¬ 
bation growth. Hence our linear calculations of erg versus 
redshift must be carefully examined and may be wrong 
due to significant non-linear corrections. 

In this vein it is interesting to note that even Qq —>■ 1 
is not strongly excluded in the condensation case. An 
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obvious test of this would be to examine galaxy rotation 
curves at low and high redshifts; one might expect the 
rotation curves at z = 0 to show the lack of dark matter. 
However, due to the nonlinear nature of the condensa¬ 
tion, the galaxies might never undergo the condensation 
process and hence might well be dark matter dominated 
despite most or all of the dark matter in the universe 
having already condensed. This nicely illustrates the po¬ 
tential subtleties involved in condensation. 

The first implication of these effects is the fact that 
metamorphosis and condensation (as we consider them) 
are likely to make the core-cusp problem of CDM worse. 

On the other hand, a recent paper [56] discussing 
Chandra observations of the elliptical galaxy NGC 4636 
finds no discrepancy between the Chandra data and the 
existence of a central dark matter cusp. The averaged 
central dark matter density inferred for this galaxy is 
about two orders of magnitude larger than for dwarf spi¬ 
rals and low-surface brightness galaxies (where the core¬ 
cusp problem seems most notable). Also intriguing, their 
fit significantly improved when using a two-component 
dark matter model rather than just a single component. 

It would also be interesting to study the formation of 
halo substructure in condensation type models. As this 
is a highly non-linear process, N-body simulations are 
needed for reliable predictions. 

A. Nonlinear effects in evaporation 

Given that the condensation and metamorphosis mod¬ 
els appear to make the core-cusp problems of standard 
ACDM worse due to the incomplete transition in regions 
of high density, it is attractive to consider evaporation 
models, as described in section (IIC). 

At the critical phase boundary the condensate, which 
exists early in the universe, starts to evaporate. At high 
temperatures, the universe consists (in the simplest mod¬ 
els) solely of baryons and the condensate in addition to 
the radiation. 

In these models, the high dark matter densities inside 
clusters would encourage a return to the condensate form 
and hence, assuming that the condensate had a long 
Compton wavelength, would suppress the formation of 
strong core density cusps which is perhaps the biggest 
current problem for ACDM. Meanwhile the IGM would 
remain in the normal, CDM, phase. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have quantitatively studied a number 
of issues in detail: 

• We have considered a phenomenological model for 
dark energy which forms via condensation of dark 
matter, which we took to be CDM. By comparing 
with current CMB, LSS and SN-Ia data we show 


that we do find a better fit to the data than from 
simple ACDM, but not significantly so. 

• The condensation model has a Compton wave¬ 
length which changes from very small to very large 
around z = Zt, in contrast to quintessence or meta¬ 
morphosis where the Compton wavelength is al¬ 
ways large. Condensation therefore has growth of 
small-scale structure absent in the other models 
which impacts on the CMB and LSS. Neverthe¬ 
less, the data slightly prefer metamorphosis, and 
strongly favour metamorphosis if Zt is forced to be 
small, Zt < 1.5. 

• We have confirmed the results of [11] for metamor¬ 
phosis even when we allow the Hubble constant to 
vary: current data favour a late-time transition in 
w(z) at Zj ~ 2 (i.e. there is no degeneracy with 
H as there is for Wf). The best-fit value of Zt for 
condensation is Zt = 4. 

We have also qualitatively discussed the nonlinear evo¬ 
lution of condensation and metamorphosis which is sig¬ 
nificantly more complex than quintessence. 

Supernova data can be used to study the evolution of 
the background geometry quite independently of the dark 
energy model. It already is able to set useful limits on 
the equation of state today. For our class of models we 
find Wf < —0.6. The supernova data also seems to prefer 
a late transition in w, but more data is needed to draw 
strong conclusions. 

Future work should consider more realistic models for 
the condensate and examine the non-linear features of 
both metamorphosis and condensation in more detail. 
Future data should improve the constraints on these 
models significantly. Eventually, the study of the proper¬ 
ties of the dark energy might lead to important insights 
into the physics at high energies which cannot be probed 
directly otherwise. 
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